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ABSTRACT 

As tracers of star formation, galaxy assembly and mass distribution, globular clusters have provided important clues to our under¬ 
standing of early-type galaxies. But their study has been mostly constrained to galaxy groups and clusters where early-type galaxies 
dominate, leaving the properties of the globular cluster systems (GCSs) of isolated ellipticals as a mostly uncharted territory. We 
present Gemini-South/GMOS g'i' observations of five isolated elliptical galaxies: NGC 3962, NGC 2865, IC 4889, NGC 2271 and 
NGC 4240. Photometry of their GCSs reveals clear color bimodality in three of them, remaining inconclusive for the other two. All 
the studied GCSs are rather poor with a mean specific frequency S'jv ~ 1.5, independently of the parent galaxy luminosity. Considering 
also previous work, it is clear that bimodality and especially the presence of a significant, even dominant, population of blue clusters 
occurs at even the most isolated systems, casting doubts on a possible accreted origin of metal-poor clusters as suggested by some 
models. Additionally, we discuss the possible existence of ultra-compact dwarfs around the isolated elliptical NGC 3962. 

Key words. Galaxies: star clusters - Galaxies: individual: NGC 3962 - Galaxies: individual: NGC 2865 - Galaxies: individual: IC 
4889 - Galaxies: individual: NGC 2271 - Galaxies: individual: NGC 4240 


1. Introduction 


The globular cluster systems (GCSs) of old, massive ellip¬ 
tical galaxies present an almost ubiquitous optical color bi - 
modality (e.g. iKundu & Whitmord l200lt iLarsen et al.l l200l1l . 
thought as corresponding to a me t allicity bimodality (e.g. 
Strader et al.l 120071: lUsher et^l2012l: iBrodie et al.ll2012[) . with 
caveats coming from possible non-linearities of the color- 
metal licity transformations introduced by horizontal-branch 


stars (lRichtledl2006l: Yoon et'^l2006tlChies-Santos et al.ll2012l: 


Blakeslee et ai.ll2oT2: Richtlerll2013ll. 

This “uniyersal” bimodality has been obyiously one of 
the main aspects the theories of GCS formation need to ad¬ 
dress. Many scenarios, not necessarily exclusiye and some- 
what oyerlapping, haye been put forward oyer the years (e.g . 
Ashman & Zenfl 1992t Forbes et al. h997^ ICote et ahl Il998 


] 11998 

Beasley e t al.l2002HRhode e t al.l2005 : lMuratoy & Gnedinj2010 : 


Elmegreen et al.ll2012l:lToninill2013). 

Since elliptical galaxies are known to inhabit mostly high- 
density enyironments (e.g. iDressled Il980ll . most of the GCS 
studies, and hence the obseryational constraints to these theories. 


* Based on observations obtained at the Gemini Observatory, which 
is operated by the Association of Universities for Research in Astron¬ 
omy, Inc., under a cooperative agreement with the NSF on behalf of the 
Gemini partnership: the National Science Foundation (United States), 
the Science and Technology Facilities Council (United Kingdom), the 
National Research Council (Canada), CONICYT (Chile), the Australian 
Research Council (Australia), Ministerio da Ciencia, Tecnologia e Ino- 
vajao (Brazil) and Ministerio de Ciencia, Tecnologia e Innovacion Pro- 
ductiva (Argentina). 


have been naturally focused on galaxy clusters le.g.lPeng et al.l 

1200^ iHarris et al.ll2006l: iMrader et al.lf2006t iLiu et al JI2()1 i1l 

In the current reigning paradigm of galaxy formation, accre¬ 
tion and merging pl ay the central roles (e.g. ICole et aLlll994t 
De Lucia et al.ll2006ll. finding strong observational suppo rt (e.g. 


Ibata et al.lll994l:lTal et al.ll2()09l: Ivan Dokkum et al.ll20roh . Sim¬ 
ulations predict a different accretion history for elliptical galax¬ 
ies in low-den sity environments compared to their high-density 
siblings (e.g. iNiemi et^ l2010ll : but the impact that a low- 
density environment may hav e on the properties of a GCS ha s 
seldom been investigated (e.g iGebhardt & Kissler-Patiglll999ll . 
although it could give imp ortant evidence t o discriminate be- 
tween accr etion driven (e.g.lCTote et al.lll99^ iHilker et"^! 19991) 
and in situ (iForbes et al.lll997 l scenarios. 

Only a handful of isolated e lliptical galaxies (lE s) have their 
GCS thoroughly studied (e.g. ISnitler et al.l l2008l) . given their 
paucity in the local environment. No f ield ellipticals are listed, 
for example, in the big compilation of iBrodie & Strad^ (120061 
see their Table 1), and only a hand ful can be seen in the newer 
compilation of iHarris et al.l (l2013l) . Recently, ICho et all d2012h 
studied 10 early-type galaxies in low density environments in the 
magnitude range -18.5 < Mg < -20 using HST/ACS, although 
solely using the lTullvl (Il988h density parameter as isolation cri¬ 
terion. 


This paper is part of an effort to understand the proper¬ 
ties of GCS in lEs and their connections to their parent galax- 
ies, collecting a large sample. This paper is a continuation of 
iLane et akl (1201 3l hereafter Paper I) who studied the GCS of the 
IE NGC 3585 and NGC 5812 using Washington photometry, and 
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Table 1. Basic data of the galaxies presented in this work. 


Name 

RA 

Dec 

1 

b 

Type 

K 

(mag) 

(m - M)q 
(mag) 

Mb 

(mag) 

fov 

(kpc) 

NGC 2271 

06:42:53.0 

-23:28:34.0 

233:13:30.5 

-12:14:45.8 

E/SO 

12.10 

32.51 

-20.41 

50.9 

NGC 2865 

09:23:30.2 

-23:09:41.0 

252:57:12.4 

+ 18:56:29.9 

E3 

12.18 

32.89 

-20.71 

60.5 

NGC 3962 

11:54:40.1 

-13:58:30.0 

282:39:10.7 

+46:38:57.5 

El 

11.59 

32.74“ 

-21.15 

56.5 

NGC 4240 

12:17:24.3 

-09:57:06.0 

289:14:11.8 

+52:00:46.3 

E/cD? 

14.01 

32.07 

-18.06 

41.5 

IC 4889 

19:45:15.1 

-54:20:39.0 

343:32:15.1 

-29:25:10.4 

E5-6/S0 

11.91 

32.33 

-20.42 

46.8 


Notes. Classification type comes from NED and HyperLeda databases, while total apparent B magnitudes are taken from RC3, except for NGC 
4240, taken from HyperLeda. Distance moduli have been adopted from iTonrv et^ l l2001h . except the one for NGC 4240 which comes from 
iReda et ^ ll2004h . The last column gives the physical size covered by the GM OS FOV at each gala xy’s distance. 

Considered as uncertain bv lTonrv et al.l (1200 ih . 


iRichtler et all ( |2015l Paper II) who studied the IE NGC 7796. It 
is als o complemented by the study of the field elliptical NGC 
7507 (ICaso et al.ll2013h also conducted by our group. 


1.1. The galaxy sample 

As isolated galaxies give the opportunity to study galaxies 
without the interference of many of the processes that af¬ 
fect their structure and evolution in more crowded environ¬ 
ments, many samples of lEs have been constructed, using 
diverse isolation criteri a ( Kuntschner et aP ItOOI ISmith et al'l 


2004IStocke et al.l2004lReda et al.l2004ICollobert et al.ll2006 


Nigoche-Netro et al.l2004 Euse et an201'^. 

In this work w e present ob s ervatio ns of five lEs tak e n from 
the catalogues of ISmith et^ (l2004h and iReda et ^ (l2004 : 
closer than ~ 50 Mpc and observable from the southern hemi¬ 
sphere. The GCS of NGC 2865 and IC 4889 have been studied 
before, albeit with a smaller field of view (see below). In the 
following we give their main characteristics as found in the lit¬ 
erature. 

- NGC 2271 is a E/S O at a surface brigh tness fluctuation dis¬ 
tance of 31.8 Mpc (ITonrv et al. 200lh presen t in the cata¬ 
logue of iReda et al.l (12004 ). Reda et akl OOOTh measured an 
age of 11.5+0.5 Gyr, with no evidence of ra dial age or metal- 
licity gradients, while iHau & Eorbei (l2006l) found strong ro¬ 
tation and kinematic signatur es of an inner disk, w hich is 
also found on optical imaging (Ide Souza et al.12004. 

- NGC 2865 is a E3 listed in the iReda et al. (l2004 cata¬ 
logue. It is known for the wealth of substructure (shells, tails, 
etc.) i t possesses dMichard & Prugnielll2004 Ide Souza et al.l 


|2004 iTal et al.l 12009 : Urrutia-V iscarra et al l |201^. being 
most probably a merger re mnant (|Hau et a D] j_999h . with a 
kinematically distinct core (iHaii & Porbesll2006h . It presents 
a remarkably steep age gradient from a centra l ~ 1.5 Gyr 
(ISanchez-Bldzauez et al.l200^lReda et al.l200^ to ~ 10 Gyr 
close to IR^; exh ibiting also a steep metal licity gradient 
( Reda et ^ l2007h . Its GCS was studied bv ISikkema et al.l 
( 2006h using HST/ACS, finding a sizeable blue population 
con sistent with a young /intermediate GC population (see 
also lTrancho et al.ll2014ii . 

NGC 3962 is part of the ISmith et^ (l2004 sample of 
held ellipticals, wi th no companions as confirmed by 
iMadore et al. ( 2004), and no signs of any tidal disturbance 
(iTal et all 20091). Annibali et al. ([2007) measured an age of 
10 Gyr, while Serra & Oosterlool ( 20101) give an age of ~2.5 
Gyr. 

NGC 4240 was taken from the [Reda et alJ (|2004 catalogue. 
iReda et ^ (l2007h found a central age of 7.5 Gyr followed 


Table 2. Observations log 


Name 

Obs. date 
(dd.mm.yyyy) 

Exp. time (s) 

g' 

EWH M (") 
g' i' 

NGC 2271 

01.02.2012 

9x419 

9x315 

0.73 

0.74 

NGC 2865 

25.01.2012 

6x480 

6x340 

0.63 

0.50 


26.01.2012 

6x490 

6x350 



NGC 3962 

22.02.2012 

6x610 

6x440 

0.73 

0.60 

NGC 4240 

02.02.2012 

5 x422 

5x215 

0.80 

0.58 

IC 4889 

17.04.2012 

3 x575 

5x355 

- 

0.75 


01.05.2012 

3x515 

1 x355 

0.49 

- 


Notes. The FWHM was measured on the combined images. 


by a mild age gradient. The existence of an also shallow 
metallicity gradient is explained as indication of a past ma¬ 
jor merger. It possesses a kinematical ly distinct core and a 
mostly flat velocity dispersion profile (iHau & Porbesll2006l) . 

- IC 4889 is a n E/SO (iLaurikainen et al.ll20111) present in the 
ISmith et ^ (|2004 catalogue, showing isophotal twisting 
(iTal et al.l 2009) as well a s a dusty disk dSmith et al.ll2004lL 
ISerra & Oosterlool (1201 Ol) estimate an age of ~ 3.6 Gyr. Its 
GCS has been studied by iGebhardt & Kissler-PatigI (1 19991) 
using HST data, without finding evidence for bimodality. 


The rest of the paper is distributed as follows; in Sect. l2.1l we 
describe the observations, as well as the data reduction and stel¬ 
lar photometry. In Sect.|3 we present the results on the color dis¬ 
tribution of globular clusters, while radial profiles are presented 
in Sect. |4] GC luminosity functions are calculated in Sect. |5] 
Sect.|6]describes properties of the galaxies themselves and their 
connections to their GCS. Sect. |7]put the results into a broader 
context, before giving the conclusions in Sect. 0 


2. Observations and data reduction 

2.1. Gemlnl-S/GMOS Imaging 

Imaging of the five galaxies was obtained using GMOS, 
mounted at the Gemini South Telescope, Cerro Pachon, Chile 
(program ID GS-2012A-Q-6). GMOS comprises three detectors 
covering in total a fleld-of-view of 5.5 x5.5 arcmin^, with a pixel 
scale of 0.146"pixel ' after a 2 x 2 binning. A set of dithered 
images in the SDSS Alters g' and i' was obtained. Details of the 
observations can be seen in Table |2] 

Image reduction was carried out inside the Gemini/iRAiQ 
package (v 1.9). Reduction included bias subtraction, flat- 

’ IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in As- 


Article number, page 2 ofM 












































































































































Salinas et al.: Globular clusters in isolated ellipticals II 


Table 3. Transformation equation coefficients for Eqs.[T]and|2]for the different galaxies. 


Galaxy 

Zero-point 
g' i' 

Color term 
g' i' 

N^stars 

rms 

g' i' 

NGC 2271 

28.226 + 0.005 

27.900 + 0.007 

-0.052 + 0.010 

-0.024 + 0.008 

11 

0.035 

0.028 

NGC 2865 

28.253 + 0.022 

27.917 + 0.012 

-0.052 + 0.010 

-0.024 + 0.008 

15 

0.046 

0.041 

NGC 3962 

28.110 + 0.034 

27.770 + 0.033 

-0.052 + 0.010 

-0.024 + 0.008 

10 

0.101 

0.103 

NGC 4240 

28.240 ± 0.008 

27.930 + 0.007 

-0.052 + 0.010 

-0.024 + 0.008 

33 

0.023 

0.021 

IC 4889 

28.209 ± 0.014 

27.923 + 0.009 

-0.007 + 0.014 

0.001 + 0.009 

20 

0.031 

0.017 


Notes. Extinction coefficients, taken from the Observatory webpages, are fixed for all equations with values Kgi = 0.18 and Kji = 0.08. Color 
terms for NGC 2271, NGC 2865 and NGC 3962, were not fitted, but fixed; see main text for details. Additionally, the last four columns indicate 
parameters for Eq.[3]for both filters as described in Sect. 12.1.21 


fielding, detector mosaicing, fringe removal (for the i' expo¬ 
sures) and image co-addition. Processed bias, flat and fringe 
frames were retrieved from the Gemini Science Archive website. 
The FWHM measured on point sources in the final images can be 
seen in Table|3 Image quality in the IC 4889 exposures showed 
considerable variations between the two dates, for this reason 
the combined g' image takes data only from the May 2012 ob¬ 
servations, while the combined /' was constructed with images 
coming from the night in April 2012 only. 

2.1.1. Stellar photometry 

In order to ease the source detection, images were first median 
filtered using a box of 10" aside. This process removes the par¬ 
ent galaxy light (everywhere but in the innermost arcseconds), 
leaving the point sources unaltered. A median value of the sky, 
determined from blank portions in the unsubtracted images dis¬ 
tant from the central galaxy, was then added back to the images 
to keep the right photon statistics during the detection and pho¬ 
tometry stages. 

Object identification was carried out using SExtractor 
(I Berlin & ArnoutsI Il996l) . The choice of the detection filter 
can have important consequences on the number of identified 
sources; while a Gaussian filter provides a better detection rate 
for the faintest sources, a “mexican hat” filter gives better re¬ 
sults for sources very close to saturated stars or to the center of 
the galaxy unsubtracted residual. Given the relatively low Galac¬ 
tic latitude of some of the targets (NGC 2271 and NGC 2865) 
and the clumpy nature of some galaxy residuals (for example, 
IC 4889, see Sect. IS, the “mexican hat” filter is the most suit¬ 
able option. SExtractor was then run twice on each image using 
both detection filters. Duplicate detections within 3 pixels were 
removed from the merged catalogues. 

Concentric aperture and psf-fitting photometr y was obtained 
using the stand-alone version of daophot/allstar (IStetsonlll987h 
using as input the object coordinates from SExtractor. Between 
30 and 50 bright isolated stars per field were chosen to deter¬ 
mine the psf shape. Aperture corrections for each frame were 
measured using a subset of these stars. 

Point sources were selected based on the behaviour of the 
SExtractor output parameter class_star and allstar param¬ 
eter sharp. The accepted range of values for both parameters 
were determined from the results of the artificial stars exper¬ 
iments (see Sect. 12.1.2b . Einal cross-search between the stel¬ 
lar positions in th e two differen t bands was done using dao- 
match/daomaster dStetsonl Il993h . keeping into the final cata¬ 
logue only sources present in both bands with a tolerance of 3 

tronomy (AURA) under cooperative agreement with the National Sci¬ 
ence Eoundation. 


Table 4. Completeness 


Galaxy 

g' 

a 

i' 

mso 

g' 

i' 

NGC 2271 

1.32 

1.63 

26.17 

25.23 

NGC 2865 

1.90 

2.95 

27.32 

26.24 

NGC 3962 

2.06 

3.04 

26.71 

25.42 

NGC 4240 

2.41 

4.24 

26.45 

25.77 

IC 4889 

3.29 

2.25 

26.82 

24.41 


pixels (~ 0.5") in their positions, thus weeding out remaining 
spurious detections. 

Standard star frames from the southern extension of the 
ISmith et^ (l2002h catalogue observed on the same nights than 
our science frames were retrieved from the Gemini Science 
Archive. These frames were reduced following the same pro¬ 
cedure of the science images, with the exception of the fringe 
correction, which was deemed to be unnecessary given the short 
exposures (~ 5 sec) in the /' filter. The number of standard stars 
measured for each night can be seen in Tabled The transforma¬ 
tion equations to the standard system used were 

9ahs = 0std +Zg + Kg(X - 1) H- bg(g' - /')std (1) 

'obs = '"td + Z/ + Ki(X-l) + biig' - i'U (2) 

where X denotes the airmass of the observation. Coefficients 
were obtained using iraf/photcal. The extinction coefficients 
were taken from the Gemini webpages and not fit, given the 
small range of airmasses covered by the observations. Also, 
when less than 20 stars were available on a field, the color terms 
were fixed to the values obtained for NGC 4240 for which the 
largest number of standard stars was available. Coefficients used 
for each galaxy can be seen in Table |3] 

2.1.2. Photometry completeness 

To study the level of completeness as a function of mag¬ 
nitude, artificial stars experiments were conducted using 
daophot/addstar. In each of the co-added frames (i.e. the ones 
where the galaxy light was not subtracted yet), we added 200 
point sources covering a wide magnitude range. These new 
frames were then processed in the same way previously de¬ 
scribed; galaxy subtraction with median filtering, source detec¬ 
tion with SExtractor and psf photometry with daophot. The 
procedure was repeated 50 times for a total of 10000 artificial 
stars per galaxy per filter. An example of the fraction of recov¬ 
ered stars as a function of magnitude can be seen in Eig. [T](top 
panel). 
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Fig. 1. Top panel: Fraction of recovered stars after the artificial stars 
experiment in the i' observations of NGC 4240. The dotted line is the 
fitted interpolation function (see Eq. [3}. The dashed blue lines indicate 
the 50% level of completeness. Middle panel: Behaviour of the all- 
star parameter sharp. Bottom panel: SExtractor stellarity index. In 
both panels artificial stars are indicated with black circles, while the 
real sources are plotted in red. The blue dashed lines indicate the point 
source criteria applied to the science data based on the artificial one. 


The fraction of recovered stars was fitted with the interpola¬ 
tion formula, 


^ 2 


1 - 


a{m - mso) 


sjl + a^(m — mso)^ 


(3) 


which gives a good analytical description of the com- 
pleteness, and have been used for many GC studies (e.g. 


McLaughlin et al.l 119941: [Fleming et al.l Il995t IWehner et al.l 
2008t Alamo-Martlnez et al.ll2012 1. The free parameters a and 


mso indicate the speed at which the detection rate falls at the 
faint level, and the 50% completeness, respectively. The values 
for both parameters for all fields and filters can be seen in Table 

m 

Artificial stars experiments are not only useful to understand 
the detection limits, but also to understand the behaviour of the 
photometric measurements at these faint magnitudes. Separa¬ 
tion between point and extended sources is a pervading prob¬ 
lem in galactic and e xtragalactic studies and several methods are 
commonly used (e.g. Harris et al.||l99l|; Bertiri & Arnoutslll99^ 


Idem et al.l200i[Durret et al.l2009HCho et al.l2012l) .'ln our case 

we preferred to use a combination of the SExtractor output pa¬ 
rameter class_star, a stellarity likelihood, and the allstar pa¬ 
rameter sharp, which measures deviations from the derived psf. 
In Fig.[Tl we show an example of the behaviour of both parame¬ 
ters as a function of magnitude for the artificial stars (black cir¬ 
cles) and the detected sources (red circles). Blue dashed lines 
show the limits that define the envelopes of the artificial stars. 
Only about l%-2% of the total artificial stars lie beyond these 
limits. In the following we consider as point sources, the sources 
that satisfy these criteria. 


2.1.3. Estimating the contamination 


GCS of elliptical g alaxies can have radial extensions up to hun - 
dreds of kpc (e.g. iRichtler et al.ll2Qlll ISchuberth et al.ll2Q12h . 
Given the relatively small FOV provided by GMOS at the galax¬ 
ies’ distance (see column 10 in Table[TJ, a control field cannot be 
take n within the same telescope pointing as i n wide-held studies 
(e.g. iRhode & ZenlbOOlUPirsch et al.ll2003L Paper I), but has to 
be selected elsewhere. Control fields are necessary in order to de¬ 
termine the contamination the GC samples suffer; contamination 
coming partly from stars and from faint background galaxies that 
cannot be culled during the point-source selection (Sect. lZTTTT i. 


Since our observing program did not include any control 
field, an appropriate field was searched in the Gemini Sci¬ 
ence Archive, finding that the William Herschel Deep Field 
(WHDF) had similar total exposure times and image quality 
than our science data, and a very similar setup, using g' and 
i' filters on Gemini- N/GMOS (Gemini ID GN-2001B-SV-104, 
iMetcalfe et al.l 120011) . The main caveat of using the WHDF as 
control field lies on its relatively high Galactic latitude (b ~ 
-61°), implying a likely underestimation of the stellar contam¬ 
ination, although at the faint levels where the bulk of the GC 
studied here are located, the dominant component of the con- 
tamin ation will be unre solved background galaxies, instead of 
stars (iFaifer et al.ll201 ih . The WHDF counts were hence com¬ 
plemented by ad ding star counts c oming from the Besan 9 on 
Galactic mo dels (iRobin et ^ 2003h at each galaxies’ position. 
As noted by iLane et al.l d201 ih . comparisons between the Be- 
san 9 on models and pencil beams surveys can be problematic, 
due to stellar streams or other Local Group objects that are not 
included in the model. In our case, neither the WHDF nor the 
studied galaxies lie close to any of these features. 


The reduction, photometry, artificial star experiments and 
point-source definition of the WHDF were canded out in the 
same way as the target galaxies, with the exception of the me¬ 
dian filtering, in the absence of a dominant galaxy. The color- 
magnitude diagram of sources in the WHDF can be seen in Fig. 
12] (bottom right panel). Since the completeness of the WHDF 
is shallower than the photometry in the NGC 2865 and NGC 
4240 fields, the contamination in the faint end of these fields 
was estimated by multiplying the number of detected contami¬ 
nant sources with the ratio of the interpolation functions (Eq. |3]l 
of the galaxies to the one of the WHDF; generating bootstrapped 
samples with the expected number of contaminant sources. 


2.2. SOAR/SAM imaging 

Additional imaging of NGC 3962 was acquired on the night of 
April 22,2014, u sing the SOAR Adapti ve Module (SAM) and its 
Imager (SAMI) (iTokovinin et al.ll201011201 3) . SAM is a ground 
layer adaptive optics system installed at the SOAR 4.1m tele¬ 
scope at Cerro Pachon, Chile. SAMI gives a field of view (FOV) 
of 3 X 3 arcmin with a scale of 0.091 "pixel ' (2x2 binning). 

In total, 5 X 120, 3 x 120 and 5 x 120 seconds exposures were 
taken in the SDSS g, r and i filters, respectively. Details of the 
reduction process will be given in a forthcoming paper (Salinas 
et al. in prep), and only a preliminary calibration to the stan¬ 
dard system, not crucial for its scientific goal, is given here. The 
FWHM measured on the combined images was 0.59", 0.55" and 
0.50" in the g, r and i filters, respectively. This set of images, of 
higher resolution than the GMOS imaging, was used to study to 
the galaxy center where the deeper GMOS images were satu¬ 
rated. More details are given in Sect. |6] 


Article number, page 4 ofIBj 

























































Salinas et al.: Globular clusters in isolated ellipticals II 


20 

21 

22 

i' 

' 023 

24 

25 


0.0 0.5 1.0 1.5 2.0 

(g'-i')o 



0.0 0.5 1.0 1.5 2.0 

(g'-i')o 





20 

21 

22 

,23 

24 

25 

26 


^ NGC 4240 • 

[Ngcs= 331 



0.0 0.5 1.0 1.5 2.0 

(g'-i')o 



Fig. 2. Color-magnitude diagrams for sources detected around the five galaxies plus the control field (bottom right). Small black dots indicate 
the resolved sources, while the big black dots are the point sources as determined in Sect. 12.1.^ Globular clusters are depicted in red, inside the 
magnitude and color limits explained in the text. The number of GC candidates stated for each galaxy is before any decontamination and within 
the FOV. Green triangles in the NGC 3962 panel indicate UCD candidates (see Section|T4). 
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3. The GC color distribution 


80 


Color-magnitude diagrams for all the point sources in our tar¬ 
get galaxies can be seen in Fig. |2] as big black points, while 
resolved sources are depicted will small black dots. Globu¬ 
lar cluster candidates were further limited to the color range 
0.4 < ig - i )o < 1.4, following previous GCS studks with the 
same filters (iForbes et al.ll2004t1wehner et al.ll2008l iFaifer et aP 
|201 Ih . but with a slightly bluer limit to a llow the study of pos - 
sible younger clusters (as in NGC 2685 ISikkema et^ 200^ 
The br ight limit was fixed at M, = -11.5 following iFaifer et al.l 
(1201 111 , and the faint limit was given by the 80% completeness 
limit in both g and i. These bounding limits can be seen in Fig. 
13 with dashed lin es. Photometry was corrected for Galactic ex- 
tinction using the Schlaflv & Finkbeinerl (1201 ll) recalibration of 
ISchlegel et al.l (11998 ) dust maps. 

The effect of the color-dependent completeness is negligible 
for NGC 3962 and the WHDF, mild for NGC 2771 and NGC 


2865 and significant for NGC 4240. There is no color-dependent 
completeness for IC 4889. 

None of the GCSs present an instantly obvious bimodal color 
distribution as t heir more massive counterparts in cluster en¬ 
vironments (e.g. iDirsch et al.ll2003h . but NGC 3962 stands out 
by its healthy GCS. The GCSs with the deepest observations, 
NGC 4240 and NGC 2865, show a significant number of blue 
ig' - i' ~ 0.5) sources below i' ~ 25, which are likely compact 
background galaxies which have survived the point source de¬ 
termination process, while IC 4889 stands out by the presence 
of a significant blue population throughout the entire magnitude 
range (21 < i < 24), which seems absent in the rest of the galax¬ 
ies (compare, for example, with NGC 3962). 

Fig. |3] shows color histograms for the five GCSs, where the 
data has been binne d using an optim um bin size, which depends 
on the sample size (llzenmannil 19911) . Each panel shows the raw 
(green dashed lines), background corrected (thick black lines) 
and background estimate (blue dotted lines) color histogram 
for each GCS. The background corrected histograms reveal that 
most of the GC candidates with colors (g' - i')o ^ 0.7 can be ac¬ 
counted as part of the contamination, except for NGC 2865 and 
especially NGC 4240 where a significant number of blue sources 
remain. T he blue populat i on in NGC 2865 is consistent with the 
results of ISikkema et al.l (l2006l) who proposed the existence of 
a population of young clusters. Visually, the clearest case for 
bimodality can be seen in the IC 4889 color histogram , which 
contradicts the result of iGebhardt & Kissler-Pati^ (Il999h which 
found no evidence of bimodality in the GCS of this galaxy. 

In order to test the existence of bimodality in the GC color 
distributions more rigorously, we use a Gaussian mixture mod¬ 
eling (gmm) as implemented bv iMuratov & GnedinI (1201 Oh . gmm 
provides three indicators useful to discriminate between uni- 
modal and bimodal distributions: the kurtosis of the distribution, 
k, which if negative indicates a flat-top distribution; the sepa¬ 
ration of the peaks, D, defined as D = l/Uwue - /^redl/[(o‘yuj. + 
which, when brger than 2, indica tes a clear sepa¬ 
ration between the peaks d Ashman et al.l[l994h : and a p-value, 
which gives the probability of obtaining the same from a 
unimodal distribution. Uncertainties in each statistic are calcu¬ 
lated with 100 bootstrapped realizations of the sample. To fur¬ 
ther test the robustness of the results, gmm was run 50 times on 
samples generated by randomly subtracting the number of con¬ 
taminants expected per each color bin. The errors measured with 
this method are always within the errors obtained from the boot¬ 
strapped realizations (which are the quoted errors in Tabled, but 
additionally provide an uncertainty for the kurtosis which is not 
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Fig. 3. Color histograms of the GCSs of the five target galaxies. The 
green dashed lines indicate the raw color distribution, the blue dotted 
line indicates the contamination in each field as determined in Sect. 
12.1.31 The thick black lines show the background-corrected distribu¬ 
tions. 
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Table 5. Globular cluster radial profiles. 


Galaxy 

power- 

O-Q 

■law fit 

n 

NGC 2271 

1.69 + 0.35 

-2.18 + 0.18 

NGC 2865 

1.38 + 0.29 

-1.88 + 0.15 

NGC 3962 

1.29 + 0.27 

-1.81 + 0.14 

NGC 4240 

-0.77 + 0.30 

-0.89 + 0.15 

IC 4889 

-1.91+0.84 

-0.61 + 0.43 


given by gmm. Given the large number of blue sources in NGC 
4240, GMM was run for sources with 0.6 < g - i < 1.4, instead of 
0.4 < g - i < 1.4 as the rest of the galaxies. 

The results of gmm applied over the backgrou nd corrected 
GC sa mples can be seen in Table |6] Following lUsher et^ 
(1201 2h we define a galaxy as showing clear bimodality when the 
conditions p-val< 0.1, Z) > 2 and k < 0 are simultaneously met. 
These criteria are clearly fulfilled by NGC 2271, NGC 2865 and 
IC 4889. NGC 4240 also fulfill the criteria, although the large 
number of blue sources makes its case less reliable. 

The most interesting case is given by NGC 3962. Despite 
being the galaxy in our sample with the richest GCS and with 
deep observations, its peak separation is barely larger than 2, 
the kurtosis is very close to zero and has the highest p-value, 
making the case for bimodality very doubtful. Intermediate- 
age globular cluster s are known to smear out bimodality (e.g. 
iRichtler et al.ll2012h . but the luminosity function shows no sign 
of intermediate-age clusters (see Sect.|5]l. Only observations us¬ 
ing a longer wavelength baseline will provide a more definitive 
answer. 

Additionally, Table|6]gives the red fraction, fr, for each GCS. 
Red fractions, simply defined as the number of red clusters over 
the total number of clusters, are further discussed in Sect. 17.21 

4. The GC radial distribution 

For photometry obtained with wide enough FOVs, the popula¬ 
tion of GCs detected around a galaxy is expected to drop off to a 
constant value at sufficiently extended radial positions from the 
galactic center. This constant number will represent foreground 
stars and faint background galaxies and the projected surface 
density of GCs as a function of galactocentric distance will give 
the radial distribution of GC candidates. In practice, given the 
rather limited FOV, we adopted the outermost bin of the radial 
distribution as background. Starting from the center of the im¬ 
ages and considering only detected GC candidates brighter than 
the 80% completeness magnitude in the i filter, we grouped GC 
candidates and contaminants into concentric circular annuli, sub¬ 
tracted the contamination in each annulus and then divided by 
the effective area over which GC candidates are distributed. 

We fitted power laws, <tqc - cror", to the surface density pro¬ 
file of galaxies in our sample, where ctgc is the number of GCs 
per square arc second and r is the projected radial distance. The 
fit was done for all GCSs in the fixed interval 40" < r < 160" 
(horizontal dashed line in Fig.lH, which avoids the central zone 
of incompleteness produced by saturation, ill-subtracted galaxy 
and possible dusty features (see Sect. Results of the fit can be 
seen as dotted blue lines in Fig. |4] and in Table |5] 

The richest GCS, NGC3962, NGC 2271 and NGC 2865 fol¬ 
low a power la w with index n -2 , as is the case for most giant 

ellipticals (e.g. lBassino et al.ll2006h . The galaxy which presents 
the larger deviation from the power law in the selected radial 
range is NGC 2865, which shows a small overabundance of GCs 







Log r [arcsec] 


Fig. 4. The radial distribution of globular clusters in the studied galax¬ 
ies. Open black triangles indicate the total population, while green filled 
squares depict the background subtracted population. Each distribution 
was fitted with a power law (blue dotted line), between the limits indi¬ 
cated by horizontal black dashed lines. 
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Table 6. GMM results. 


Galaxy 

N 

/^blue 

7* red 

<^blue 

<^red 

fr 

D 

K 

p-val 

NGC 2271 

295 

0.896 ±0.014 

1.184 ±0.025 

0.125 ±0.011 

0.067 ±0.011 

0.17 ±0.06 

2.87 ± 0.26 

-0.43 ± 0.04 

0.001 

NGC 2865 

332 

0.848 ± 0.043 

1.125 ±0.044 

0.128 ±0.019 

0.096 ± 0.020 

0.38 ±0.14 

2.43 ±0.19 

-0.87 ± 0.03 

0.001 

NGC 3962 

557 

0.947 ± 0.057 

1.208 ±0.041 

0.156 ±0.019 

0.089 ± 0.019 

0.19 ±0.12 

2.05 ± 0.37 

-0.11 ±0.04 

0.007 

NGC 4240 

123 

0.765 ± 0.055 

1.098 ±0.063 

0.108 ±0.033 

0.103 ±0.033 

0.44 ±0.17 

3.08 ±0.35 

-1.16 ±0.04 

0.001 

IC 4889 

74 

0.922 ±0.015 

1.187 ±0.019 

0.062 ± 0.010 

0.081 ±0.010 

0.54 ± 0.07 

3.65 ± 0.43 

-1.25 ±0.05 

0.001 


Notes. The second column indicates the number of (background subtracted) GC candidates used as input for GMM. Columns 3 and 4 indicate the 
peak color of the blue and red subpopulations, while columns 5 and 6 indicate the dispersion of the best fitting Gaussians. Column 7 indicates the 
red fraction, that is, the number of cluster associated with the red subpopulation divided by the total number of GC candidates. Column 8 is the 
peak separation as defined in Sect. [3] while column 9 is the sample kurtosis. 


Table 7. GCLF fitted parameters and specific frequencies. 


Galaxy 


m 

O-m 

Ngc 

Sn 

Tn 

Tblue 

Led 

NGC 2271 

84 ± 1 

24.54(f) 

1.16 ±0.02 

562 ±9 

1.46 ±0.38 

5.87 ±0.10 

4.87 ± 0.25 

1.00 ±0.35 

NGC 2865 

61 ±2 

24.92 (f) 

1.02 ± 1.1 

410 ± 8 

1.00 ±0.19 

2.51 ±0.05 

1.55 ±0.48 

0.95 ± 0.35 

NGC 3962 

128 ± 15 

24.50 ± 0.21 

1.44 ±0.15 

854 ± 98 

1.26 ±0.57 

2.87 ± 0.32 

2.33 ± 0.39 

0.55 ± 0.35 

NGC 4240 

- 

- 

- 

84 ±9 

2.15 ±0.80 

2.53 ± 0.27 

1.42 ±0.55 

1.11 ±0.45 

IC 4889 

- 

- 

- 

280 ± 17 

0.85 ±0.12 

2.55 ±0.15 

1.17±0.19 

1.38 ±0.20 


at log r ~ 1.95 (90") which might be associated with its intricate 
shell system. 

NGC 4240 and IC 4889 present very shallow GC profiles, 
probably an effect of the low number of GCs present in these 
systems. Particularly problematic is the case of IC 4889 where 
an abrupt jump in the cluster distribution is visible at r ~ 40". 
IC 4889 not only possesses a poor cluster system, but also is the 
galaxy with the shallowest observations of the set. 

Finally, we note that the galaxies with fewer dusty features 
(NGC 2271 and NGC 4240, see Section O are the ones which 
follow closer the power law distribution inside the fit limit of 
40". 


5. Luminosity functions and distances 


Luminosity functions in the i' band for the point sources in 
each galaxy field can be seen in Fig. |5] (blue dotted lines). 
Background/foreground contamination was estimated using the 
WHDF and Besan 9 on models as in Sect. 12.1.31 (green dashes 
lines). The background-corrected luminosity functions, that is, 
the ones for the globular cluster candidates, are depicted with 
solid blue lines. 

The inverted red arrows indicate the expected position of 
the turnover magnitud (TOM) of the globular cluster luminos¬ 
ity fun ctions (GCLFs), assuming a universal value of = 

-8.46 ( Kundu & Whitmor^l2001 ). transformed to M, using Eq. 
2 from iFaifer et al.l(l2011() . and distance moduli as given in Table 
□ 


For NGC 2271, NGC 2865 and NGC 3962, the 
contamination-subtracted GCLFs were modelled as the 
product of a Gaussian function with the completeness function 

{Eqm- 


dN flo 

(m — nio)^ 

2 

^ a(m - mso) 


2cr^ 

L ^ m J 


^Jl + a^(m - m^o)^ 


(4) 

where mo is the TOM and cr„, the width of the GCLF. The pa¬ 
rameters for the completeness interpolation function are the ones 
derived in Sect. 12.1.21 The fitting of equation |4] was done us- 
ing MPFiT, an impl ementation of non-linear square fitting in IDL 
dMarkwardtl l2009l) . The luminosity functions were fitted down 


to the 50% completeness detection limit in /' (Table |4|i, with 
the exception of NGC 2865 were the fit was done using magni¬ 
tudes down to the 80% limit, given the large number of sources 
at i > 25.5, which are probably faint background galaxies that 
escaped the point-source criteria defined in Sect. 12.11 

NGC 3962 is the only galaxy where the completeness limit 
is well beyond the TOM, therefore only for this galaxy we fit mo, 
cTm and the Gaussian amplitude, uq. In the cases of NGC 2271 
and NGC 2865, the TOM was fixed t o its expected value based 
on its SBF distance dTonrv et al.l200lh . and only cr^ and ao were 
fit. The results of the fit for these three galaxies can be seen in 
Table 12] 


Only the case of NGC 3962 provides material to discuss its 
distance based on its GCLF. The distance m odulus derived from 
the TOM position is 32.47 + 0.20, while the lTonrv et al.l d2001h 
value, assumed throughout this paper, and based on SBF mea¬ 
surements is 32.74 + 0.40. Even though they agree within the er- 
ros, the GCLE based distance is significantly lower than the SBE 
distance. If the SBE distance were the correct one, this would 
imply a slightly brighter TOM, TOM, = -8.24, instead of the 
expected -7.97. 


As discussed in Sect. [3] a possible explanation for the ab¬ 
sence of a clear separation between the color peaks in NGC 3962 
is the presence of an intermediate age population, but intermedi¬ 
ate age populations have exactly the opposite effect, producing 
fainter TOMs (e.g. lRichtlerll2003h . also even though metallicity 
is expected to have an influence on t he TOM, this i s minimized 
when using I band observations (e.g. Reikuba[2012h . A brighter 
TOM agrees with iBlakeslee & TonrvI (1 19961) who found evi¬ 
dence that the TOM varies with environm ent, with low-densit y 
galaxies having a brighter TOM (see also IVillegas et al.ll2010l) . 
A larger sample of lEs with deep enough observations would be 
necessary to confirm this. 


The radial profiles in Eig. |4|indicate that the density of GCs 
quickly drops to values close to zero within the EOV, therefore, 
we did not attempt an extrapolation of the total number of clus¬ 
ters, Ncc, to obtain the fraction outside the EOV. The total num¬ 
ber of GCs is then simply derived from the Gaussian fitted to 
the GCLE. In the cases where no fitting was attempted (NGC 
4240 and IC4889), the total number of GCs was estimated by 
doubling the number of detected clusters above the TOM. The 
GCLE in NGC 4240 differs from the expected shape close to a 
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Fig. 5. GCLF for the five GCSs. Blue lines indicate the raw GCLF 
(dotted), the LF of the contamination (dashed), and the contamination- 
corrected GCLF (solid). Red lines indicate the fitted function (Eq.|4l 
dotted line) and its respective Gaussian. Vertical dashed black lines in¬ 
dicate the 80% and 50% completeness limits in i' as obtained in Sect. 

[2T21 


Gaussian, with a large number of very faint sources. Given the 
old age of the galaxy (Sect. Ol and the lack of evidence from 
the color maps of a recent meger event (Sect.|6]), these are prob¬ 
ably a residual contamination from background galaxies and not 
a faint young populatio n of clusters. A s imilar GCLF is seen in 
NGC 3377 as shown bv ICho et al.l(l2012h . 

Total number of GCs for the five galaxies can be seen in 
TableE] 


6. Galaxies surface brightness and colors 

Surface brightness profiles of the galaxies were measured via 
ellipse fitting using the iraf task ellipse. All the other objects 
in each galaxy field were initially masked using the segmen¬ 
tation image produced by SExtractor, but additional regions 
(extended halos of bright stars, chip defects, etc) were interac¬ 
tively masked during the ellipse execution. Once the regions to 
mask were established, ellipse was run a second time, fixing the 
center of the ellipses as the mean central value of the ellipses 
with semi-major axes between 30 and 300 pixels. Sky bright¬ 
ness was calculated using the unmasked p ixels outside the l ast 
measured isophote, taking the robust mean (iBeers et al.lll990l) of 
40 random samples containing 800 sky pixels each. The uncer¬ 
tainty in the sky value was taken as the dispersion of these 40 sky 
measurements. Photometric calibration was done with the same 
equations as the point source ph otometry (Table [S, correcting 
for Galactic extinction following ISchlaflv & Finkbeineil (l201 II) 
dust maps. 

The results from the ellipse fitting can be seen in Appendix 
IS The errors in the g' and /' surface brightness profiles comes 
from the sum in quadrature of the rms scatter in the intensity 
from the ellipse fitting and the uncertainty on the sky value. Sur¬ 
face brightness is given as function of the geometric mean ra¬ 
dius, r - a Vl - e, where a is the fitted ellipse major axis and e 
its ellipticity. 

Color maps were also generated for each galaxy in order to 
look for color anomalies related to dust and possible past merg¬ 
ers. 

NGC 2271 shows a flat color profile in the inner ~ 30" along 
the major axis, with noticeable color gradient at larger radius 
(Fig IXB. This is the norm al behaviour for early-type galaxies 
(e.g. E a Barbera et alJl201^ . and com patible with t he neg ligible 
age and metalicity gradients found by Reda et ^ (l2007n inside 
Q.SRg. No dusty features are visible. 

Fig. IA.2I (left panel), shows the color map of NGC 2865. 
NGC 2865 is widely considered as a merger remnant (iHau et al.l 
Il999l) . The right panel shows the galaxy surface brightness map 
with the same scale. The large of number of s hells, bluer than 
the galaxy body, have been discussed elsewhere (iMalin & Cartel! 
Il983l : l^rt et alJ[T986l) . We just emphasize a very narrow ~ 20" 
long and ~ 1.5" wide str eam (white box in Fig. lA.21 i. within the 
structure labeled “jet” in iFort et~^ dl986h (their Fig. 3), point¬ 
ing in projection directly to the galaxy center. Its blue color to¬ 
gether with the large shell at the exact opposite side of the galaxy, 
probably makes it part of the same accretion process that pro¬ 
duced the shells, although these very narrow radial features are 
not reproduced b y simulations (e.g. iHernauist & OuinnI [l987l : 
iBIlek et ^l2014t) . which show rather fan-like structures. 

Fig. lA.3l shows th e q - i color map fo r NGC 3962. The same 
structures detected bv iBuson et al.l(Il993h using Ha imaging can 
be seen. First, an extended, slightly off-centered arc-like struc¬ 
ture about 15" from the galaxy center, with redder color than 
the underlying galaxy light. The SAM images (Fig. IA.31 right 
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panel), also show a central elongated component, misaligned 
with the major and minor axis of the galaxy. Ionized gas is now 
commonly thought to be of external origin , especially for mas¬ 
sive, round ellipticals (e.g. lSarzi et al.ll2()()^ . so even though the 
large scale structure of the galaxy does not reveal any evidence 
for past accr etion, as indicat ed by the very low tidal parameter 
measured bv lTal et^ (l2009h . the inner color structure might be 
indication of an older merger. 

NGC 4240 shows a flatter color profile (Fig. lA.4ll. similar to 
the on e measured on the field elliptical NGC 7507 ( Salinas et al.l 
I 2 OI 2 I) . T his agrees w i th the shallow age and metallicity gradients 
found bv iReda et al.l (l2007h . A notable feature is the A{g - i) ~ 
0.15 color difference, for the inner ~ 15", between the South and 
North sections. Even though one woul d be tempted to relate it to 
the lop-sided rotation curve found by iHau & ForbesI (l2006l) on 
the same scales, that rotation pattern was found roughly along 
the East-West axis. Another feature is a red plume of ~ 15" in 
the SE direction (indicated with a black arrow in Fig lA.41 i. Both 
features were undetected bv iReda et al.l (l2004h . 

IC 4889 shows a large spiral dusty feature (Fig. lA3T l. very 
similar to the one on NGC 3962. A remarkable feature is the 
unphysically red external halo. Despite our efforts to understand 
any issues with these images, we were unable to pinpoint the 
origin for the very flat i profile. 


7. Discussion 


7 .1. The specific frequency of dusters in iow- and 
high-density environments 

The globular cluster specif ic frequency, S n, defined by 
Harris & van den Berghi(ll981l) as 

SN ^ (5) 

connects the luminosity of a galaxy with its globular cluster pop¬ 
ulation, in an attempt to measure how efficient is globular cluster 
formation compared to star formation. In Eq. |5] Ncc is the total 
number of globular clusters in a galaxy, and My is its absolute V 
magnitude. 

We obtain S n values for our sample using Ncc as obtained 
from the GCLF analysis in Section |5] Figure |6] (upper panel) 
contains a comp arison between our results (red stars) with the 
ICho et all (1201 2h sample of 10 low-density early-type galaxies 
(blue squares) and galaxies from the ACS Virgo Clus ter Sur¬ 
vey (ACSVCS, open black triangles. iPeng et Sl2008h . Addi¬ 
tionally, we include the of five more isolated/lo w-density el- 
lipticals, NGC 35 85, NGC 5 812 dLane et al.ll2013h . NGC 7507 
(ICaso etalj2013h. N GC 720 (iKartha et al.ll20l4 and NGC 7796 
( Richtler et alJl201^ (green filled circles). 

Fig. |6] (upper left panel) presents 5 at as a functio n of stellar 
mass. Stellar masses for each galaxy in our sample and ICho et alJ 
(1201 2h sample were derived using homogeneous Ks apparent 
magnitudes taken from the 2MASS Extended Source Catalog 
(iJarrett et al.ll2003l) . /T-band luminosities were transformed into 
stellar mass es by using th e Hubb le type depe ndent mass-to-light 
ratios from ISnitler et al.l (l2008h . based on iBruzual & Chariot 
(l2003h population models. Stellar masses for the ACSVCS 
galaxies were taken directly from IPeng et al.l (l2008l) . Distances 
for the ACSVCS galaxies come from the reca librated surface 
bright ness fluctuations measurements given by iBlakeslee et al.l 

Isolated ellipticals have a low Sn, scattered around Sn ~ 1-5, 
although there is a trend in the 10.5 < log M* < 11.4 in the sense 
that less massive galaxies have S n closer to 2, while the more 



Fig. 6. Specific frequency (upper panels) and globul ar clus t er red frac¬ 
tion (lower panels) as function of stellar mass and the lTullvl l ll988h den¬ 
sity parameter. Re d symbols r eprese nt our sample, while blue squares 
are taken from the ICho et all i2012f) sample. Bla ck open triangles are 
the results from the ACSVCS dPmg et aill2008h . while the the green 
filled circles are five more IBs described in the text. 


massive are clustered around ~ 1. A Spearman’s rank correlation 
test gives a coefficient of -0.74, with a two-sided significance of 
0.002, indicating a high correlation. Considering a similar mass 
range, the o pposite trend is s een for the Virgo GCSs (already 
discussed bv iPeng et alJl2008h . lEs with logM* >11.4 show a 
diverse behavior, while NGC 3962 follo ws the trend of less mas- 
sive galaxies, NGC 720 (S n = 3 + 0.2 . iKartha et al.ll201^ and 
NGC 7796 {Sn = 2.5 + 0.5. iRichtler et alj|2015h show an in¬ 
creased S M, although still lower than Virgo ellipticals of similar 
mass. 

The upper right p anel of Fig. |6] shows now the 5 ^ as a func¬ 
tion of iTullvl (119881) density parameter. Galaxies in both low- 
density and high density environments share a broad range of 
S AT, but again a trend within the low-density sample can be seen, 
where the most isolated galaxies (p ~ 0.1) have even lower S m, 
altough this trend almost totally driven by NGC 7507, the galaxy 
with the lowest density in the sample and with a surprisingly 
low Sn (ICaso et al.N2013h . Note that several galaxies in Virgo 
lack a density parameter, as well as NGC 4240 and NGC 7796, 
which are not included in the iTullvl (Il988h catalog due to their 
distances. 

The S ff depends on My, which in turns depends strongly on 
the age. Given that isolated ellipticals ar e considered late accre- 
tors, spanning a w ide range of ages (e.g. lKuntschner et al.ll2002l : 
[Reda et ^l2007h . a true comparison with cluster ellipticals, be¬ 
lieved to be homogeneously old, is hampered. To compare galax¬ 
ies with possible different stella r populations, it is releva nt to use 
the T AT parameter introduced bv IZepf & AshmanI (Il993h . 


Tn = Nccl(MJl0\ (6) 

where the number of GCs is normalized to the stellar mass, M*, 
instead of My. Additionally, the parameters Tbiue and T^ed can 
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Fig. 7. T -parameters (see text) as function of stellar mass and density 
parameter. Symbols are the same as in Fig|^ 



be defined, where the number of blue and red clusters are used 
instead of the total number, Ncc- The amount of blue and red 
clusters are obtained using the red fractions measured in Sect. [3] 

Figure |7] shows Tbiue and Fred for our sample tog ether with 
the sa me galaxies used in the S ^ comparison; the ICho et al] 
(1201 2h sample of low-density environment ea rly-type galaxies, 
and early-type galaxies from the ACSVCS (iPeng et al.l l2006i 
l2008h . Colors and symbols are the same as in Fig.|6] 

The trend with mass seen for the Sff is not seen for Tbiue, 
although the more massiye isolated ellipticals still show signif¬ 
icantly lower yalues when compared to Virgo ellipticals of the 
same mass. The Tred yalues for our sample are systematically 
lower than the ICho et all ( 1201 2b and ACSVCS samples in the 
range of masses described. 

Finally, the right panels of Fig. |7] shows the behayior of the 
Tied and Tbiue parameters as function of the density parameter. 
The trend seen for S n among isolated and low-density galaxies 
is repeated for Tbiue although less clearly. No clear trend is seen 
for Treri, although if we solely considering our sample without 
the ICho et all (1201 2b galaxies is clear than the T^ed yalues for 
isolated ellitpicals span a narrow range close to 1 or less, while 
Virgo ellipticals haye on ayerage higher yalues. 


7.2. The red fractions of isolated ellipticals 

The proportion between blue and red clusters i n a GCS is usually 
measured using the red fraction of GCs (e.g. lPeng et ahll^Obb . 
that is, the fraction of red clusters compared to the entire popula¬ 
tion in a GCS. Red fractions are an important tool to discriminate 
between the possible scenarios of GCS formation. 

Red fractions for our galaxies deriyed in Sect.[3]are presented 
in Fig. |6] (lower panels, red stars), tog ether with red fra ctions 
for GCSs belonging to the ACSVCS (iPeng et al.l 120081 black 
triang les) and the 1 0 low-density ell ipticals from ( Cho et aP 
1201 2l blue squares). iPeng et ^ (l2008b red fractions ar e an up - 
dated yersion of the results presented by iPeng et aTI (l2006b . 
where red fractions were corrected for the radial incomplete¬ 


ness that affects mostly the higher-luminosity ellipticals in their 
sample. Uncertain t ies ha ye been taken from the original re- 
sults of Peng et al] (l2006b since updated yalues are not giyen by 
|Peng_et^ ( 2 008b. U ncertainties for the red fractions measured 
by Cho et al.l (1201 2b are also not giyen in their paper. 

The red fractions of Vi rgo galaxies prese nt the same non¬ 
monotonic behayior seen in iPeng et al.l (l2008i their Figure 8) as 
function of instead of stellar mass. Red fractions increase 
from lower mass galaxies up to a maximum around log M* ~ 11 
beyond which red fractions experience a turnoyer. Interestingly, 
the somew hat lower red fr actions for low-density ellipticals 
claimed by lCho et al. ( 2012b a re not seen; in the 10 < logM* < 
11 range the ICho et al ( 2012b sample is indistinguishable from 
the Virgo g alaxies. That we cannot see lower red fractions for the 


ICho et alT(l2012b samp le compared to the ACSVCS is probably 
based on, a) ICho et"^ (1201 2b did not use the red fractions cor¬ 
rected by aperture giy en by|Peng et ah ( 2008b . instead using the 
uncorrected ones from IPeng et al.l (l2006b . and b) the comparison 
of red fractions was done as a function of Mb instead of stellar 
mass; for the same stellar mass, low-density ellipticals will ap¬ 
pear brighter in the optical as they are on ayerage younger than 
cluster ellipticals, so the claimed lower red fractions are only a 
population effect on the galaxies luminosities and not an intrinsic 
difference between the GCSs. 


Howeyer, the galaxies presented in this paper appear to show 
systematic lower fractions than the Virgo galaxies at similar stel¬ 
lar mass, with the exception of IC4889. IC 4889 has the shallow¬ 
est obseryations in our sample, thus the red fraction obtained for 
its GCS would appear to be uncertain, accounting for its unusu¬ 
ally high yalue. It is important to note that, just like the yalues for 
the total number of GCs, the fr are not corrected for unaccounted 
clusters outside the FOV, based on the quick decline of the GC 
density profiles (Sect. nil. If anything, wider field studies of these 
isolated galaxies will reyeal eyen lower red fractions. We there¬ 
fore conclude that red fractions in isolated ellipticals are equal 
or possibly lower than the yalues for high-density ellipticals. 

Recently, iToninil (1201 3b introduced a model of GC forma¬ 
tion based on hierarchical merging and the mass-metallicity re¬ 
lation for galaxies. An important feature of the model is the abil¬ 
ity to make predictions for galaxies with rich and poor accretion 
histories, where poor-accretion galaxies end up with dominant 
red GC populations. If we associate galaxies with poor accre¬ 
tion histories to present day isolated ellipticals, our results are 
in stro ng disagreement wit h this model and in line with the re¬ 
sults of iKartha et al.l(l2014b for lenticular galaxies. Additionally, 
it would be inconsistent with models proposing an a ccreted ori¬ 
gin fr om dwarf galaxies for the blue population (e.g. lCdte et al.l 
Il998h . 

A strong argument against the accreted origin of blue clus¬ 
ters is the correlation between the peak color of the blue clusters 
and the luminosity of the parent galaxy (e.g. iLarsen et al.ll200ll: 
IStrader et al.l l2004b . which is not expected under an accretion 
scenario. Our hnding of a dominant blue peak despite the likely 
absence of a rich accretion history, fayours in-situ formation for 
GCs, where then the color bimodality would be rathe r explained 
by a non-lin ear color-metallicity relation in GCs (e.g. lYoon et al.l 
1200611201 lb . 


7 . 3 . The color peaks 

The mean color, or color peak, of the blue and red GC subpopu¬ 
lations haye bot h been found to be functions of the paren t galaxy 
luminosity (e.g. iLarsen et al.ll200ll IStrader et al.ll2004b . which 
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Fig. 8. Black stars indicate the red and blue peaks of the col or distribu¬ 
tion of GCs in our sample. Squares indicate galaxies from the lCho et alJ 
( 12012|) sample, w hile open trinagles are the values from the ACSVCS 
dPeng et al.ll2008h . 


suggest an intimate connection between galaxy and GC enrich¬ 
ment history. _ _ 

In order to compare again with ICho et af] (1201 2h and 
iPeng et aTI(l2008ll samples, we transform (g-z) color s from those 
studie s into (g - i) using the relation obtained by lUsher et^ 
(1201 2h for 169 bonafide cluster in NGC 4365, 


(g-i) = (0.735 + 0.009) x(g-z)H-(0.147+ 0.012) mag. (7) 

Fig. [8] shows that on average both peak magnitudes in our 
sample ar e redde r, b y a small but signi ficant amount, than the 
ICho et all (1201 2h and IPeng et al.l ( ^08h peak m agnitudes. This 
is in contrast to the results of Cho et al] (1201 2h who found that 
the peak magnitudes for t heir GCS were slig htly bluer compared 
to the Virgo galaxies from IPeng et al.l (l2008l) . which also fits into 
the picture of isolated/field ellipticals having in average younger 
stellar populations. Given that our redder values go against this 
intuition, we believe our data may be affected by an error in its 
zeropoint photometry. We therefore restrict the analysis to an 
internal comparison between the galaxies in our sample. 

Even though for the blue clusters, a slope may be more no¬ 
ticeable, this is dominated by the most massive (NGC 3962), 
and least massive galaxies (NGC 4240) of our sample, which 
are precisely the ones for which color bimodality is not totally 
supported (see Sect. O. Moreover, the mass-peak color relation 
for blue clusters is exp ected to be even les s pronounced for low- 
density environments (IStrader et al.ll2004h . We conclude that no 
relation with galaxy mass is discernible for our sample, mostly 
limited by the narrow range of masses studied. 


7.4. Ultra-compact dwarfs in NGC 3962? 

The CMD of point sources around NGC 3962 shows a clear 
extension of the GC sequence above the selected brighter limit 
(Fig. |2] top left panel, green triangles). Seventeen point sources 
brighter than i - 21.24 (M, = -11.5 at the galaxy’s distance) 
and with 0.8 < (g - i)o <1.4 can be seen. In this same range 
of luminosity and colors only about three point sources are ex¬ 
pected based on the comparison field (the WHDF, Fig.|2] bottom 


right plot) and the Galactic model of iRobin et al.l (l2003h . Stellar 
systems of this bightness hav e been commonly labele d as “ultra¬ 
compact dwarfs” (UCDs, e.g. jDrinkwater et^l2000h . due to the 
absence of counterparts in the Local Group. 

Even though some of the low luminosity UCD candidates 
can be just bright GCs, given the uncertainty of 0. 4 mag in the 
surfac e brightness distance modulus measured by iTonrv et aH 
(1200 Ih . the correspondence between the expected and measured 
peak of the GCLF in NGC 3962 (Fig. |5]) shows that the true dis¬ 
tance cannot be too off, and hence most of these are genuinely 
brighter than the expected brightest GCs. 

UCDs have been found mostly in galaxy gro ups and clus¬ 
ters (e.g. iMieske et al.|l2008t iBrodie et aTllj^Ollh . with forma¬ 
tion channels mostly appealing to dynamical effects common 
in high-density environments (e.g. iHilker et al.Nl999f) . If spec¬ 
troscopically confirmed, the UCDs in NGC 3962 would be the 
first around an isolated elliptical, confirming that the formation 
of these objects is not restricted to high-density environments, 
followi ng the finding of one UCD around the field spiral galaxy 
M 104 (!Hau et al.ll2009h . 

8. Summary and conclusions 

In this paper we have studied the GCSs of five isolated elliptical 
galaxies, NGC 2271, NGC 2865, NGC 3962, NGC 4240 and IC 
4889, using Gemini-S/GMOS gi imaging. Our main results are: 

- The GMM analysis indicates clear bimodality for three out 
of 5 galaxies, while the two remaining show very short peak 
distances, making their bimodality less clear. In general, 
color bimodality is a common phenomenon down to the most 
isolated elliptical galaxies. 

- The specific frequency of isolated ellipticals, S n, is close 
to 1.5 irrespectibly of the galaxy mass. In the range 10.5 < 
log M* < 11.4 a small correlation can be seen, in the sense 
that more massive isolated ellpiticals have smaller S n, but 
for higher masses, larger S m a re me asured. 

- The conclusion of ICho et al.l (1201 2h that galaxy mass is the 
main ruler of the properties of a GCS stems from the nar¬ 
row range of masses used in their sample, biased towards 
low-luminosity galaxies. For galaxies with stellar masses 
log M* >11, environment, and not mass, has the highest im¬ 
pact determining some of the global properties (e.g. specific 
frequencies) of the GCSs. 

- The red fractions show no significant difference with the 
ones seen in ellipticals in Virgo, making an accreted origin 
for the bl ue clusters doubtf ul. The slightly lower red fractions 
found by ICho et all (l2Q12h are a stellar population effect on 
the galaxies luminosities, and not intrinsic to the GCSs. 

The last point assumes isolated ellipticals can be connected 
to poor accretion history galaxies found in cosmological sim¬ 
ulations. Even though many isolated ellipticals show signs of 
recent mergers, to acquire knowledge of the entire accretion his¬ 
tory, that is, to establish whether present day isolated ellipticals 
have been is olated throughout the ir lifetimes, is more difficult. 
As shown by ICooner et al.l (1201 3l) , the outer surface brightness 
profiles of high-density and low-density galaxies should be re¬ 
markably different given their different accretion histories. The 
study of the outer surface brightness profiles of isolated ellipti¬ 
cals as proxy of their accretion histories will be the subject of 
forthcoming publications. 
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Fig. A.l. NGC 2271 color map. 



Fig. A.4. NGC 4240 color map. 



Fig. A.5. IC 4889 color map. 


Appendix B: The surface brightness profiies 

g and i surface brightness profiles are given for each galaxy as 
measured with iraf/ellipse. The indicated radii are the mean geo¬ 
metric radii of each fitted ellipse. Radii are given in arc seconds, 
while surface brightness, in mag/arcsec^. 


Appendix A: The coior maps 

g - i color maps are given for all five galaxies. Each panel has 
a different angular scale, appropriate to indicate relevant color 
features in each galaxy. For all figures North is up. East, to the 
left. 
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Fig. A.2. Left panel'.g — i color map of NGC 2865. Right panel: g surface brightness map of NGC 2865. Both figures share the same size and scale. 
The black/white box indicates the radial feature discussed in the text. 
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Fig. A.3. g - i color map for the central parts of NGC 3962. Left panel: Gemini/GMOS map, 1' aside. The central hlue area is due to detector 
saturation. Right panel: SOAR/SAM color map, 30" aside. Albeit with lower 5/A, it shows details unseen in the Gemini image due to saturation. 
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Table B.l. NGC 2271 surface photometry 


Table B.2. NGC 2865 surface photometry 


R 

9 

?err 

i 

^err 

R 

9 

^err 

/ 

^err 

2.36 

18.630 

0.000 

17.352 

0.000 

3.65 

18.557 

0.000 

17.547 

0.000 

2.58 

18.719 

0.000 

17.441 

0.000 

3.99 

18.681 

0.000 

17.665 

0.000 

2.82 

18.810 

0.000 

17.531 

0.000 

4.37 

18.808 

0.000 

17.784 

0.000 

3.08 

18.902 

0.000 

17.621 

0.000 

4.77 

18.938 

0.000 

17.906 

0.000 

3.36 

18.995 

0.000 

17.715 

0.001 

5.22 

19.073 

0.000 

18.038 

0.000 

3.66 

19.091 

0.000 

17.813 

0.001 

5.72 

19.210 

0.000 

18.172 

0.000 

4.00 

19.189 

0.000 

17.912 

0.001 

6.20 

19.340 

0.000 

18.305 

0.000 

4.35 

19.287 

0.000 

18.013 

0.001 

6.77 

19.476 

0.000 

18.444 

0.000 

4.73 

19.387 

0.000 

18.114 

0.001 

7.40 

19.621 

0.000 

18.581 

0.000 

5.16 

19.495 

0.000 

18.223 

0.001 

8.09 

19.764 

0.000 

18.727 

0.000 

5.63 

19.612 

0.000 

18.341 

0.001 

8.85 

19.905 

0.000 

18.871 

0.000 

6.16 

19.738 

0.000 

18.469 

0.001 

9.65 

20.039 

0.001 

19.001 

0.000 

6.74 

19.869 

0.000 

18.597 

0.001 

10.53 

20.168 

0.001 

19.136 

0.000 

7.35 

19.997 

0.000 

18.731 

0.001 

11.56 

20.308 

0.001 

19.274 

0.001 

8.03 

20.131 

0.000 

18.869 

0.001 

12.68 

20.458 

0.001 

19.424 

0.001 

8.79 

20.273 

0.000 

19.012 

0.002 

14.02 

20.633 

0.001 

19.594 

0.001 

9.58 

20.413 

0.000 

19.156 

0.002 

15.52 

20.819 

0.001 

19.785 

0.001 

10.48 

20.569 

0.001 

19.315 

0.002 

17.17 

21.011 

0.001 

19.987 

0.001 

11.46 

20.733 

0.001 

19.482 

0.003 

18.95 

21.212 

0.002 

20.189 

0.001 

12.55 

20.911 

0.001 

19.664 

0.003 

21.00 

21.428 

0.002 

20.410 

0.001 

13.77 

21.103 

0.001 

19.861 

0.004 

23.19 

21.648 

0.002 

20.632 

0.002 

15.08 

21.295 

0.001 

20.054 

0.004 

25.91 

21.888 

0.003 

20.874 

0.002 

16.66 

21.522 

0.001 

20.288 

0.005 

28.99 

22.126 

0.004 

21.121 

0.003 

18.44 

21.747 

0.002 

20.521 

0.007 

31.63 

22.310 

0.004 

21.310 

0.003 

20.44 

21.976 

0.002 

20.760 

0.008 

34.64 

22.519 

0.005 

21.522 

0.004 

22.93 

22.226 

0.002 

21.013 

0.011 

38.28 

22.767 

0.006 

21.762 

0.005 

25.75 

22.483 

0.003 

21.288 

0.014 

42.17 

22.987 

0.008 

21.983 

0.006 

28.98 

22.732 

0.004 

21.550 

0.017 

46.38 

23.164 

0.009 

22.175 

0.007 

32.47 

22.969 

0.005 

21.795 

0.022 

51.02 

23.363 

0.011 

22.378 

0.009 

36.39 

23.211 

0.006 

22.058 

0.028 

56.12 

23.554 

0.013 

22.583 

0.011 

40.43 

23.435 

0.007 

22.305 

0.035 

61.73 

23.758 

0.016 

22.793 

0.013 

44.47 

23.649 

0.009 

22.520 

0.042 

67.91 

24.007 

0.020 

23.026 

0.016 

48.92 

23.868 

0.011 

22.747 

0.052 

74.70 

24.300 

0.026 

23.333 

0.021 

53.81 

24.080 

0.013 

22.988 

0.064 

82.17 

24.616 

0.034 

23.671 

0.029 

59.19 

24.316 

0.016 

23.223 

0.079 

90.39 

25.054 

0.051 

24.099 

0.043 

65.11 

24.572 

0.021 

23.500 

0.101 

99.42 

25.436 

0.072 

24.586 

0.067 

71.63 

24.781 

0.025 

23.663 

0.117 

109.37 

25.627 

0.086 

24.885 

0.087 

78.79 

24.898 

0.028 

23.885 

0.142 

120.30 

25.880 

0.107 

25.081 

0.104 

86.67 

25.201 

0.036 

24.123 

0.174 

132.33 

26.171 

0.140 

25.564 

0.159 

95.33 

25.349 

0.042 

24.288 

0.200 






104.87 

25.547 

0.050 

24.375 

0.215 






115.35 

25.558 

0.050 

24.390 

0.218 
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Table B.3. NGC 3962 surface photometry 


Table B.4. NGC 4240 surface photometry 


R 

9 

?err 

i 

^err 

R 

9 

fcr 

i 

^err 

5.82 

18.990 

0.000 

17.844 

0.001 

1.23 

17.672 

0.000 

17.039 

0.000 

6.38 

19.116 

0.000 

17.976 

0.001 

1.34 

17.780 

0.000 

17.042 

0.000 

6.98 

19.239 

0.000 

18.104 

0.001 

1.47 

17.915 

0.000 

17.045 

0.000 

7.63 

19.361 

0.000 

18.226 

0.001 

1.61 

18.051 

0.000 

17.075 

0.000 

8.34 

19.478 

0.001 

18.344 

0.001 

1.76 

18.183 

0.000 

17.180 

0.000 

9.13 

19.594 

0.001 

18.469 

0.001 

1.92 

18.313 

0.000 

17.305 

0.000 

10.10 

19.725 

0.001 

18.595 

0.001 

2.09 

18.438 

0.000 

17.427 

0.000 

11.10 

19.858 

0.001 

18.740 

0.001 

2.26 

18.554 

0.000 

17.542 

0.000 

12.29 

20.010 

0.001 

18.900 

0.002 

2.44 

18.664 

0.000 

17.650 

0.000 

13.61 

20.177 

0.001 

19.074 

0.002 

2.63 

18.766 

0.000 

17.753 

0.000 

15.07 

20.355 

0.001 

19.257 

0.002 

2.84 

18.863 

0.000 

17.858 

0.000 

16.59 

20.530 

0.001 

19.440 

0.002 

3.08 

18.971 

0.000 

17.972 

0.000 

18.26 

20.710 

0.002 

19.630 

0.003 

3.35 

19.091 

0.000 

18.102 

0.000 

20.06 

20.891 

0.002 

19.808 

0.003 

3.71 

19.235 

0.000 

18.256 

0.000 

21.97 

21.061 

0.002 

19.989 

0.004 

4.21 

19.417 

0.000 

18.449 

0.000 

24.17 

21.239 

0.003 

20.176 

0.005 

4.90 

19.631 

0.000 

18.656 

0.000 

26.66 

21.430 

0.003 

20.366 

0.006 

5.54 

19.835 

0.000 

18.852 

0.000 

29.33 

21.612 

0.004 

20.553 

0.007 

6.15 

20.006 

0.000 

19.031 

0.000 

32.23 

21.787 

0.004 

20.725 

0.008 

6.79 

20.170 

0.000 

19.202 

0.000 

35.10 

21.936 

0.005 

20.869 

0.009 

7.53 

20.343 

0.000 

19.380 

0.000 

38.18 

22.069 

0.006 

21.000 

0.010 

8.32 

20.515 

0.000 

19.555 

0.000 

41.88 

22.218 

0.006 

21.150 

0.012 

9.13 

20.689 

0.000 

19.732 

0.000 

46.26 

22.375 

0.007 

21.319 

0.014 

10.02 

20.893 

0.000 

19.937 

0.000 

51.01 

22.547 

0.009 

21.499 

0.016 

11.03 

21.117 

0.000 

20.167 

0.000 

55.94 

22.723 

0.010 

21.678 

0.019 

12.16 

21.344 

0.000 

20.399 

0.000 

60.98 

22.882 

0.012 

21.841 

0.022 

13.32 

21.546 

0.000 

20.603 

0.000 

66.48 

23.067 

0.014 

22.031 

0.027 

14.62 

21.736 

0.000 

20.798 

0.000 

73.13 

23.282 

0.017 

22.257 

0.033 

16.04 

21.917 

0.000 

20.985 

0.001 

80.44 

23.526 

0.021 

22.489 

0.040 

17.58 

22.092 

0.000 

21.174 

0.001 

88.49 

23.775 

0.027 

22.744 

0.051 

19.19 

22.243 

0.001 

21.311 

0.001 

97.33 

24.028 

0.033 

22.921 

0.059 

20.99 

22.415 

0.001 

21.523 

0.001 

107.07 

24.378 

0.046 

23.360 

0.088 

23.13 

22.602 

0.001 

21.688 

0.001 

117.78 

24.712 

0.062 

23.611 

0.109 

25.34 

22.780 

0.001 

21.880 

0.001 

129.55 

25.094 

0.086 

23.771 

0.125 

27.97 

22.996 

0.001 

22.091 

0.001 

142.51 

25.323 

0.107 

24.316 

0.201 

30.94 

23.256 

0.001 

22.384 

0.002 

156.76 

26.018 

0.196 

25.150 

0.396 

33.94 

23.521 

0.002 

22.650 

0.002 

172.43 

26.870 

0.395 

26.267 

0.875 

37.39 

23.813 

0.002 

22.960 

0.003 






41.13 

24.108 

0.003 

23.220 

0.004 






45.50 

24.411 

0.004 

23.618 

0.006 






50.50 

24.732 

0.005 

23.930 

0.008 






55.89 

25.041 

0.007 

24.277 

0.011 






61.95 

25.388 

0.009 

24.630 

0.015 






69.30 

25.777 

0.013 

24.969 

0.021 






74.52 

26.024 

0.017 

25.146 

0.025 
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Table B.5. IC 4889 surface photometry 


R 

9 

?err 

i 

^err 

2.57 

17.980 

0.000 

17.319 

0.000 

2.83 

18.092 

0.000 

17.319 

0.000 

3.12 

18.203 

0.000 

17.319 

0.000 

3.41 

18.324 

0.000 

17.320 

0.000 

3.74 

18.453 

0.000 

17.372 

0.000 

4.11 

18.581 

0.000 

17.502 

0.000 

4.52 

18.716 

0.000 

17.635 

0.000 

4.98 

18.849 

0.000 

17.768 

0.000 

5.46 

18.980 

0.000 

17.904 

0.001 

6.00 

19.123 

0.000 

18.042 

0.001 

6.67 

19.298 

0.000 

18.209 

0.001 

7.28 

19.428 

0.000 

18.338 

0.001 

8.03 

19.574 

0.000 

18.482 

0.001 

8.75 

19.699 

0.000 

18.677 

0.001 

9.65 

19.837 

0.000 

18.830 

0.001 

10.64 

20.009 

0.000 

18.930 

0.001 

11.70 

20.160 

0.000 

19.075 

0.001 

12.79 

20.303 

0.001 

19.216 

0.002 

13.96 

20.444 

0.001 

19.354 

0.002 

15.27 

20.594 

0.001 

19.502 

0.002 

16.78 

20.764 

0.001 

19.651 

0.003 

18.32 

20.927 

0.001 

19.812 

0.003 

19.91 

21.083 

0.001 

19.967 

0.003 

21.75 

21.247 

0.001 

20.135 

0.004 

23.87 

21.437 

0.002 

20.316 

0.005 

26.52 

21.664 

0.002 

20.508 

0.006 

29.27 

21.891 

0.002 

20.718 

0.007 

32.38 

22.132 

0.003 

20.932 

0.008 

36.16 

22.403 

0.004 

21.168 

0.010 

40.36 

22.687 

0.005 

21.411 

0.013 

44.99 

22.965 

0.006 

21.638 

0.016 

50.01 

23.222 

0.008 

21.845 

0.019 

55.68 

23.493 

0.010 

22.005 

0.022 

61.74 

23.748 

0.013 

22.235 

0.027 

68.51 

24.027 

0.017 

22.463 

0.033 

79.26 

24.486 

0.025 

22.616 

0.038 

88.12 

24.755 

0.032 

22.798 

0.045 

97.77 

25.153 

0.046 

23.011 

0.054 

110.10 

25.633 

0.071 

23.285 

0.069 

122.16 

26.178 

0.115 

23.613 

0.093 

139.19 

27.188 

0.269 

23.959 

0.126 
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